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1. Introduction 

The correlation functions of the spin- 1/2 XXZ spin chain have attracted 
much interest in mathematical physics through the last two decades. One 
5J] ' of the most fundamental results is the exact derivation of their multiple- 

5t 1 integral representations. The multiple-integral representations of the XXZ 

correlation functions were derived for the first time by making use of the 
g- vertex operators through the affine quantum-group symmetry in the mas- 
sive regime for the infinite lattice at zero temperature.^''^ They were also 
derived in the massless regime by solving the q-KZ equations. ^'^ Making 
use of algebraic Bethe-ansatz techniques such as scalar products, ^"^" the 
multiple-integral representations were derived for the XXZ correlation func- 
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tions under a non-zero magnetic field. "'^^ They were extended into those at 
finite temperatures/^ and even for a large finite chain. ^"^ Interestingly, they 
are factorized in terms of single integrals. ^^ Furthermore, the asymptotic 
expansion of a correlation function of the XXZ model has been systemat- 
ically discussed. ^^ Thus, the exact study of the XXZ correlation functions 
should play an important role not only in the mathematical physics of in- 
tegrable models but also in many areas of theoretical physics. 

The Hamiltonian of the spin- 1/2 XXZ spin chain under the periodic 
boundary conditions is given by 



1 ^ 
Hxxz - ^ ^ {^f<^f+i + '^J'^J+i + Acrf crf+i) • (LI) 



i=i 

Here cr" (a = X, Y, Z) are the Pauli matrices defined on the jth site and A 
denotes the XXZ coupling. We define parameter q by 

A = iq + q-')/2. (1.2) 

We define rj hy q = expvy. In the massive regime: A > 1, we put rj = 
( with C > 0. At A = 1 (i.e. q — 1) the Hamiltonian (1.1) gives the 
antiferromagnetic Heisenberg (XXX) chain. In the massless regime: — 1 < 
A < 1, we set ?7 = i(, and we have A — cosC with < C < ^i" for the 
spin- 1/2 XXZ spin chain (1.1). In the paper we consider a massless region: 
< C < 7r/2s for the ground-state of the integrable spin-s XXZ spin chain. 

Recently, the correlation functions and form factors of the integrable 
higher-spin XXX and XXZ spin chains have been derived by the algebraic 
Bethe-ansatz method. ^'^®"^® The solvable higher-spin generalizations of the 
XXX and XXZ spin chains have been derived by the fusion method in 
several references. ^^"^''' In the region: < C < t^/'^s, the spin-s ground-state 
should be given by a set of string solutions. ^®'^^ Furthermore, the critical 
behavior should be given by the SU(2) WZWN model of level k = 2s with 
central charge c = 3s/(s + i). 24, 28-40 p^j. ^j^^ integrable higher-spin XXZ 
spin chain correlation functions have been discussed in the massive regime 
by the method of g- vertex operators. "^^^^^ 

In the present paper we discuss several important points in the algebraic 
Bethc-ansatz derivation of the correlation functions for the integrable spin-s 
XXZ spin chain where s is an arbitrary integer or a half-integer.^ In par- 
ticular, we briefly discuss a rigorous derivation of the flnite-sum expression 
of correlation functions for the spin-s XXZ spin chain. 

The content of the paper consists of the following. In section 2 we formu- 
late the i?-matrices in the homogeneous and principal gradings, respectively. 
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They are related to each other by a similarity transformation. In section 3 
we introduce the Hermitian elementary matrices and construct conjugate 
basis vectors for the spin-s Hilbert space in the massless regime. In section 
4 we construct fusion monodromy matrices. In section 5, we first present 
formulas^ for expressing the Hermitian elementary matrices in terms of 
global operators. Then, we review the multiple-integral representations of 
the spin-s XXZ correlation function for an arbitrary product of elementary 
matrices.^ In section 6 we briefly sketch the derivation of the finite-sum 
expression of correlation functions for the spin-s XXZ spin chain, which 
leads to the multiple-integral representation in the thermodynamic limit. 
Here the spin-1/2 case corresponds to cq. (5.6) of Rcf. 11. 

2. Symmetric and asymmetric i?-matrices 

2.1. R-matrix and the monodromy matrix of type (1, 1®^) 

Let us now define the i?-matrix of the XXZ spin chain. ^"^'^^ For two- 
dimensional vector spaces Vi and V2, we define R^{Xi — X2) acting on 

Vi (g) V2 by 

/I 0\ 



i?±(Ai-A2)= y i?±(u)^^e"'^®e^''^ 



,.fc,c,d=0.1 



b{u) c=f(m) 
Oc±(m) b{u) 

Vo 1/ 



(2.1) 
where u = Ai — A2, b{u) = sinhM/sinh(u -I- 77) and c^{u) = 
cxp(±u) sinhry/ sinh(M + rj). We denote by e°' ^ a unit matrix that has only 
one nonzero element equal to 1 at entry (a, b) where a, 6 = 0, 1. 

The asymmetric i?-matrix (2.1), R^lu), is compatible with the homo- 
geneous grading of Uq{sl2)-^'^^ We denote by R'^p^u) or simply by R{u) 
the symmetric i?-matrix where c (u) of (2.1) are replaced by c{u) = 
sinhr// sinh(M + ?7).^* It is compatible with the affine quantum group Uq{sl2) 
of the principal grading. ^'^^ Hereafter, we denote them concisely by R^'^\u) 
with w = ± and p, where w — + and w = p in superscript show the homo- 
geneous and the principal grading, respectively 

Let s be an integer or a half-integer. We shall mainly consider the tensor 
product Vi (^ ■ ■ ■ (^ Vj^f ^ of (2s + l)-dimensional vector spaces V " 
with parameters ^j, where L = 2sNs. Here Ng denotes the lattice size of 
the spin-s chain. In general, we may consider the tensor product Vjy ^^ (E) 
Vi " (g) • • • (X) Iv ^'' with 2si + • • • + 2sr = L, where V * have parameters 
Xj or ^j for j = l,2,...,r. For a given set of matrix elements A^^'^ for 
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a, 6 = 0, 1, ... , 2sj and a, /3 = 0, 1, . . . , 2sk, we define operator Aj^k by 

a,fc=l a,^ 

(2.2) 

Here E'^' ^^ denote the elementary matrices in the spin-Sj representation, 
each of which has nonzero matrix element only at entry (a, b). 

When So ~ £/2 and si ~ ■ ■ ■ Sr — s, we denote the type by (i, (2s)®^=). 
In particular, for s ~ 1/2, wc denote it by {£, 1®^). 

2.2. Gauge transformations 

Let us introduce operators $j with arbitrary parameters 0^ for j = 
0, 1, ... ,L as follows: 

^^-(IX) -^«'^'^«(L">^"''"^'^- (2.3) 

In terms of Xjk = '^j^ki we define a similarity transformation on the K- 
matrix by 

^ffe = X,fei?,fcX7fc' ■ (2.4) 

Explicitly, the following two matrix elements are transformed. 

(^?.)'^ = c(A„ A,)e*^-^'= , (i?^,)^' = c(A„ A,)e-^^+^'= . (2.5) 

Putting 0j = Aj for j = 0, 1, . . . , i in cq. (2.3) we have 

^%{^3.^^) = iXjkf^ R,k{^j.\k) {XjkV^ (j,fc = 0,l,...,i). (2.6) 

Thus, the asymmetric i?-matrices Ri2{^1t^2) are derived from the sym- 
metric one through the gauge transformation Xjk- 

2.3. Monodromy matrices 

Applying definition (2.2) for matrix elements Riu)"^^ of a given _R-matrix, 
i?("''(u) for w == ± andp, we define i?-matrices R^ {\j, Afe) = Rjk i-^j^^k) 
for integers j and k with < j < fe < L. For integers j, k and £ with 
Q<j<k<£<L, the i?-matrices satisfy the Yang-Baxter equations 

^jk i-^3 ^ ^k)Rj^ (Aj - Xi)Rm i^k - ^e) 
= ^u (^fe ~ ^i)Rj^ (Aj - ^e)Rj^ (Aj - Afe) . (2.7) 
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Let US introduce notation for expressing products of i?-matrices. 

^1,23-n — -n-ln ' ' ' ^13 -"-12 ' 
-n-12-n-l,n — -n-ln -n-2n '''^n-ln- \^-°) 

Here R^ denote the _R-matrix R^' = R^' {Xa — Xb) for a,b ^ 1,2, ... ,n. 
We now define the monodromy matrix of type (1, l^^w), i.e. of type 
(1, 1"^^) with grading w. Expressing the symbol (1, 1®^) briefly as (1, 1) in 
superscript we define it by 

^0,^12-^1 (-^o; {wj}l) = RqlO'O - Wl) • • • Ro2i>'0 - W2)i?cJi(Ao - Wi) 
= ^OL ^OL^l ■ ■ ■ ^01 ^ Ro^l2-Li^o;{Wj}L)- (2.9) 

Here we have put A^ = Wj for j = 1,2, ... ,L. They are arbitrary. We call 
them inhomogeneous parameters. We express the operator-valued matrix 
elements of the monodromy matrix as follows. 

°'''-'' ' \cil']^{X-,{w,},)D^^+},{X-,{w,w) 

We also denote the operator-valued matrix el- 

ements by [Tq ij-M'^! {^i}i)]a-b f°'' o,-,^ = 0,1. Here {wj}l denotes the 
inhomogeneous parameters wi,W2, . . . ,wl- Hereafter we denote by {jjlj^m 
the set of N numbers or parameters jii, . . . , iij^. 

The monodromy matrix of principal grading, Tq x2-''-l('^' {^j}i)i i^ '"®" 
lated to that of homogeneous grading via similarity transformation xoi-i = 
$0^1 • • • ^L as follows.^® 

^0,\2.^l(A; {w^j}l) = X012...LTo*,\2.^.i(A; {w-j}l)Xo12-L 

^ ( Xi2-LA[\''l^iX;{w,}L)Xi2'..L e-^"xi2-LBg''UX;{w,}L)xl2-L\ 
[e^»Xi2-LC[l''}^{X;{w,}L)x^l.^ Xi2-lD[1'UX;{w,}l)x^I,^ J 

In Ref.^* operator A(^+'(A) has been written as j4+(A). 

2.4. Operator R: Another form of the R-matrix 

Let Vi and V2 be (2s + l)-dimensional vector spaces. We define permutation 
operator Hi_ 2 by 

Hi, 2 vi iSi V2 = V2 iSi vi , wi e Fi , W2 e V2 • (2-11) 

In the spin-1/2 case, we define operator .rJ™^j^(u) by 

i?(;;ViH = %.+i4!^>iN- (2-12) 
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3. The quantum group invariance 

3.1. Quantum group Uq{sl-2) 

The quantum algebra Uq{sl2) is an associative algebra over C generated by 
X^^K^ with the following relations:'*^"'*'^ 

[X\X-] = —--^. (3.1) 

The algebra Uq{sl2) is also a Hopf algebra over C with comultiplication 

A(X+) =X+ ®l + K®X+ , A{X-) =X- (g) R-^ + 1(^X- , 
A{K)^K®K, (3.2) 

and antipode: SiK) = K'\S{X+) = -K-^X+ ,S{X-) = -X-R, and 
coproduct: e(X±) = and e{K) = 1. 

It is easy to see that the asymmetric i?-matrix gives an intertwiner of 
the spin-1/2 representation of Uq{sl2): 

i?]*2(w)A(x) ==To A(x)i?+(M) for x = X^,K. (3.3) 

Here we remark that spectral parameter u is arbitrary and independent of 
X± or K. 



3.2. Temperley-Lieb algebra 

Operators Rj jj^iiu) are decomposed in terms of the generators of the 
Temperley-Lieb algebra as follows.' 

R%+i{^)=I-Kn)Uf. (3.4) 

U^s {U~s) satisfy the defining relations of the Temperley-Lieb algebra:'^ 

ufuf_,,uf = uf, 

Uf+i^f^t+i = Uf^ f"^-?' = 0, 1, . . . , i - 2, 

{Uff = {q + q-^) uf for.?- = 0, 1, . . . , L - 1, 

Ufuf - Ufuf for \j - fc| > 1 . (3.5) 

We remark that the asymmetric i?-matrices R, ,_|_i(u) derived from the 
symmetric i?-matrix through the gauge transformation are related to the 
Jones polynomial.''^ 
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3.3. Basis vectors of spin-ill representation of Uq{sl2) 

Let us introduce the g-integer for an integer n by [n]q = {q"' — q^"') / (q—q^^)- 
We define the q- factorial [n]ql for integers n by 

[n]q\ = [n]q [n - 1], • • • [1], . (3.6) 

For integers m and n satisfying m > n > wc define the g-binomial coeffi- 
cients as follows 

We now define the basis vectors of the {£ + l)-dinicnsional irreducible 
representation of Uq(sl2), \\f^,n) for n — 0, !,...,£ as follows. We define 
11^,0) by 

11^0) = 10)1® |0)2®---®|0),. (3.8) 

Here \a)j for a = 0, 1 denote the basis vectors of the spin-1/2 representation 
defined on the jth position in the tensor product. We define | \£, n) for n > 1 
and evaluate them as follows^* . 

IIA") = (a(^-i)(x-))"|K,o)^ 

Y^ cr-...^-|0)g»i+»2+-+»"-«^+«("-i)/2. (3.9) 
l<n<-<i„<£ 
We define the conjugate vectors explicitly by the following: 

'- -'9 l<ii<---<in<l 

(3.10) 
It is easy to show the normalization conditions:^® (^j'^ll IK:*^) = 1- In the 
massive regime where q = exprj with real 77, conjugate vectors (^,fi|| are 
Hermitian conjugate to vectors ||^, n). 

3.4. Conjugate vectors 

In order to construct Hermitian elementary matrices in the massless regime 
where \q\ = 1, we now introduce another set of dual basis vectors. For a 
given nonzero integer i we define (^, f^| | for n — 0,1, ...,n, by 

M=(^) E (0|a+...a+<7-(*i+-+^")+"^-"("-i)/2. 

^ ^ l<n<-<i„<£ 

(3.11) 



\n\gl 
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They are conjugate to ||^, n): (i?,m|| ||^, n) = 5m,n ■ Here we have denoted 
the binomial coefficients for integers £ and n with < n < ^ as follows. 

(3.12) 



nj [i — n)\n\ 

We now introduce vectors | \i.^ n) which are Hermitian conjugate to (^, n| | 
when |g| = 1 for positive integers d with n = 0, 1, . . . ,i?. Setting the norm 

of ||£,n) such that (£,n|| ||£,n) = 1, vectors \\t,n) are given by 

y^ err ...crr|0)g"(*i+---+*"'+"-^""(""i)/2 ^ g-»(^-«)/^^^ 

^ ...^ ^ ^ C LJ/7 \/ 



l<ii<-<i„<£ 



<? 



(3.13) 
We have the following normalization condition: 



(^,n||||£,n) = 



(3.14) 
n I 

g ' 



3.5. Hermitian elementary matrices 

In the massless regime we define elementary matrices E™' " '-^^ +■* for m,n = 
0,l,...,2sby 

j^m, „ (2. +) ^ 1 1-^^^ ^2s, n| I . (3.15) 

In the massless regime where \q\ = 1, matrix ||^, n)(^, n|| is Hermitian: 
(11^, n)(^, n||)'l' = ||£, n)(i?,n||. However, in order to define projection opera- 
tors P such that PP = P, we have formulated vectors \\i, n). 

3.6. Projection operators 

We define projection operators acting on the 1st to the fth tensor-product 
spaces by 

n'2^.-Ell^'")(^'"ll- (3.16) 

n=0 

Let us now introduce another set of projection operators -P|...^ as follows. 

£ 

Pfl, = Y.\^){i,n\\. (3.17) 

n=0 



May 7, 2010 0:51 WSPC - Proceedings Trim Size: 9in x 6in DM-RVF 



Projector P^.. ^ is idempotent: (P;{...^)^ — P\..i- In the massless regime 

where \q\ = 1, it is Hermitian: (i\...^) — Pi...£- From (3.16) and (3.17), 
we show the foUowing properties: 

pW pW _ pW ^o 1 s,\ 

^12--l^l---e. ~ ^12-£ ' l^O.iOj 

P\^.P\S... = P^.- (3.19) 

In the tensor product of quantum spaces, V^ ® ■ ■ ■ ® V^ * , we define 

mL by 

P^■L = \{PU-^n^■ (3.20) 

i=l 

Here we recaU L = 2sNs- 

The projection operators are also constructed by the fusion method. For 
i > 1 we can construct projection operators inductively with respect to t 
as follows.25'46 

Pfl, = nt.''-i4ti,,((^ - l)^)At.'.^ ■ (3-21) 

The projection operator P12...1 gives a q-analogue of the full symmetrizer 
of the Young operators for the Hecke algebra.'*® 



4. Fusion construction 

4.1. Higher-spin monodromy matrix of type {I, {2s)®^') 

We now set the inhomogeneous parameters Wj for j ~ 1,2, . . . , L, as Ng sets 
of complete 2s-strings.^* We define Wfj-iU+fl for /3 = 1, . . . , 2s, as follows. 

"'SLi)+0=^b-(/3-l)^- for 6- 1,2,..., TV,. (4.1) 

We shall define the monodromy matrix of type (1, (2s)®^=) associated with 
homogeneous grading. We first define the massless monodromy matrix by 

'f(l,2s+) Q . r^ 1 N _ p(2s) p(l.l+)('\ ./„,, (2^5)1. ^p(2s) 

^ /I(2-^+)(A;{6}ivj|(^*+HA;{6}ivj\ .,.. 
l^C(2^+)(A;{6}A'Ji?(2^+HA;{6}A^Jy' '^ ■ ^ 

Let us introduce a set of 2s-strings with small deviations from the set 
of complete 2s-strings. 



and/3== l,2,...,2s. (4.3) 



<s(b-i)+0-^b-(/5-l)^ + eCS for 6 = l,2,...,iV, 
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Here e is very small and r^ are generic parameters. We express the el- 
ements of the monodromy matrix T^^'^' with inhomogeneous parameters 
given by w\ *' "^ for j = 1, 2, . . . , i as follows. 



T, 



(1,1+), 



0, 12---L 



„(2«;^)\,\- 



A 



(2s 
12---L 



>iX)B\ 



(2s 



[\{Wj iL)- „(2s+;£).,. ri(^s+,e) 



'(A) 



a 



12---L 



'{X) D^X'W 



(4.4) 



Here A 



(2s+;< 
12---L 



'(A)denotes42:.^^(A;{wP^^}L 



(2s; e)^ 



A 



'(2s+) 
12---N. 



(A;feVJ = lin.F(^14l^:l^)(A;{^f -) 



}l)pI 



(2s) 
2---L ■ 



(4.5) 



We define the massless monodromy matrix of type {£, {2s)'^'^=) by 

n(l,2s+) ., N^(l,2s+) 
ai, 12---N^y'^ai)-l-a2, 12---N 

X • • • T^:'l2.^1r, (Aa, -{£- 1)V) P^'\,a.-a, ■ (4.6) 






4.2. Integrable spins Hamiltonians 

We define the massless transfer matrix^ of type {£, {2s)^'^') by 

?2."^t^ (A) - tr^,, (f^^^^ (A)) = E « (^' -II ri:; ^X W X 

n=0 
X T^i:;?2.tk(^ - ^) • • • ^Sl2.S.(A - (^ - 1)^) K^)a ■ (4-7) 

It follows from the Yang-Baxter equations that the higher-spin transfer 



(2s), 



matrices commute in the tensor product space V\ 
derived by applying projection operator Pi2^..i^ to V\ 



V^ * , which is 



V, 



(1) 



The massless spin-s i?-matrix R^ 



(2s, 2s +) 



(u) becomes the permutation 



operator at u = 0: R12 * (0) = ni^2-^^'^° Therefore, putting inhomo- 
geneous parameters ^p = for p = 1, 2, . . . , A^^, we show that the trans- 
fer matrix 42^'w (A) becomes the shift operator at A = 0. We derive 
the massless spin-s XXZ Hamiltonian by the logarithmic derivative of the 
massless spin-s transfer matrix. 



-l/(2s) _ « , 7(2s,2s+).,^ 



A=0,{j=0 






(4.8) 



M=0 
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5. Spin-£/2 massless XXZ correlation functions 

5.1. Spin-s local operators in terms of global operators 

In the massless regime, we can express the Hermitian elementary matrices 
in terms of global operators as follows.^ For m > n we have 



E] 



m,n (1+) 



^r-' 



pi^.L n (^^'+'+^^'+^)K) 



fe— 1 /c— n+l A;— m+1 

X n (A(i+)+i?(i+))K)pW^. (5.1) 

For m < n we have 

(i-i) 



s: 



■m, n (^ +) 



I 



fc=i 



fc— m+1 



fe=n+l 



X n (yl(i+)+i?(i+))K)p(!'.^. (5.2) 

By the quantum inverse-scattering problem (QISP) of Rcf. 10 the local 
spin operators are expressed in terms of global operators and the transfer 
matrices for the integrable spin-s XXX spin chain. However, it is not clear 
how one can derive (5.1) and (5.2) by the QISP method even for q = l. 

5.2. Symbols for expressing sequences 

Let us denote by {aj)m a sequence of numbers aj for j = 1,2, ... ,m, i.e. 
(aj)m — (ai, ^2, • • • ,0.m)- 

Definition 5.1. We say that a sequence (6fe)„ is a subsequence of (aj)m 
if (i) n < m, (ii) b^ € {ai, . . . , Cm} for k = 1, 2, . . . , n, (iii) for any pair of 
integers j and k satisfying 1 < j < fc < n, there exists a pair of integers 
^(j) and £(fc) such that aj = ^^(j), cLk — &£(*;) and (i{i) < £{k). 

For a pair of sequences {aj)m and {hk)n, wc define the product 
{o-j)m^{bk)n by a sequence {c^)m+n such that Cj = aj for j = 1, 2, • • • ,to 
and Cj = bj for j = m + l,m + 2, 



,m + n. 
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5.3. Conjecture of the spins Ground-state solution 

Let us now introduce the conjecture that the ground state of the spin-s 
case \tpg " ) is given by Ns/2 sets of 2s-strings: 

A)i") ^ fia-ia-l/2)rj + ei°'\ for a= 1,2,..., Ns/2 and a == 1,2,..., 2s. 

(5.3) 
Here we assume that string deviations ei" are very smaU when Ng is very 
large. '^ In terms of Xa , the spin-s ground state in the homogeneous grading 
is given by-'^ 

Ns/2 2s 

l^f ^^) = n Y{B^^'^H\^a^-{S.v}Nm. (5.4) 

a— 1 Q — 1 

We denote by M the number of Bcthc roots: M = 2s Ns/2 ~ sNg. 

According to analytic and numerical studies, ^'*''^^''^^'*° we may assume 
the following properties of string deviations ea s. For very large Ns, the 
deviations are given by ea = i Sa , where i denotes \/—l and Sa are real. 
Moreover, 6^^^ - <5i"+'^ > for a ^ 1, 2, . . . , 2s - 1, and \6t^\ > |,5i"+'^| 
for a < s, while \Sa \ < \5a \ for a> s. 

In the limit: Ng — > oo, the density of string centers, ptotifJ-), is given by 

1 ^' 1 

Ptot{^J■) = T7 y" :^7 rr-f rTTTT ■ (5-5) 

^'f^i 2Ccosh(7r(^-Cp)/C) 

For the homogeneous chain where ^p = for p = 1,2, . . . , Ng, we denote 
the density of string centers by p(A) . 

P^^^ = 2Ccosh(7rA/C) • (^-^^ 

Let us introduce useful notation of the suffix of rapidities. For ra- 
pidities Xa = A(a,a) wc define integers ^ by A = 2s(a — 1) + a for 
a — 1,2, . . . ,Ns/2 and for a — 1,2,..., 2s. We thus denote X(^a,a) also 
by A^ for ^ = 1,2,..., sNs, and put X(^a.a) in increasing order with respect 
to A = 2s{a- 1) +a such as A(i_i) = Ai, A(i^2) = A2, . . . , A(jv,/2,2s) = A^jv,- 
In the ground state, rapidities Xa for A = 1, 2, ... , M, are expressed by 

A2s(a-i)+a = Ma - (a - 1/2)7? + £<") (1 < a < Ng/2; l<a< 2s). (5.7) 

For A = 2s(a — 1) + a with 1 < a < 2s, integer a is given by a = [{A — 
l)/2s] + 1, and integer a is given hy a = A — 2s[{A — l)/2s]. 

For a real number x we define [x] by the greatest integer less than or 
equal to x. We define a{j) and a{j) for j = 1, 2, . . . , M as follows. 

a{j)^[{j-l)/2s] + l, a{j)=j-2s[{j^l)/2s]. (5.8) 
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5.4. Correlation functions of the integrable spin-s XXZ 
model on a long finite chain 

We define the correlation function of the integrable spin-2s XXZ spin chain 
for a given product of (2s + 1) x (2s + 1) elementary matrices such as 

^n,ji ( s+) . . .^^'J™ ( '*~''^ on the spin-s ground state, {ipg "^ ), as follows. 

fc=i 

(5.9) 
By formulas (5.1) and (5.2) we express the mth product of (2s + 1) x (2s + 
1) elementary matrices in terms of a 2s?Tith product of 2 x 2 elementary 
matrices with entries {cj, £,} as follows. 

m 2sm , 

I[Er'''"'^'-Ci{^,MPtA • I[er^ . P^!,. (5.10) 

6=1 k=l 

By making use of (5.1) and (5.2), C{{ib,jb}) is given by 

I . (5.11) 



C({^fc,Jfc}) = ^ (^' 

6=1 l^-^'' 



2s 
ib 



2s 
jb 



Here e2s(6-i)+/3 and £2^(6-1)+^ (6 = 1, . . . , iV^; /? = 1, • • • , 2s) are given by 

_/! (l</3<i6) ' _(l (1</S<*6) 

£2.(6-1)+/^ - \ (76 < /3 < 2s) ' 'Mb-i)+f3 - I {lb < 13 < 2s). 

(5.12) 
We evaluate the spin-2s XXZ correlation function F,„ ({**;, Jfc}) by 

2sm f 

X n^"^ • /^l2!'j^f+V(^f +^l^f +^>(5.13) 
i=i 

Let a"'' be the set of j with ej — 0, and a^ the set of j with e. = 1: 

a+-{j;e,=0}, a" = {j; e; = 1} . (5.14) 

We denote by r and r the number of elements of the set ex and a+, 
respectively. Due to charge conservation, we have 

r + r = 2sm. (5.15) 

We denote by jmin and jmax the smallest element and the largest element 
of a~, respectively. We also denote by j^-^^ and jmax the smallest element 
and the largest element of a^, respectively 
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Recall that the ground state \ipg ) has M Bcthe roots with M — sNg. 
Let Cj (j E a^) and c- {j E «+) be integers such that 1 < Cj < M for 
j E a^ and 1 < c- < M + j for j E «+. We define sequence (bi)2sm by 

(61, 62, . . . , b2sm) = i(^j'^^^^ ■ • • ' S"' - ' '^J"'i"' • • • ' '^J""-) • (^--^^^ 



Here sequence (c ., , . . . , c ., , c^^;^ , • • • , Cj„^^ ) is given by the composite 

sequence of c- s in decreasing order with respect to suffix j, and Cj s in in- 
creasing order with respect to suffix j. We introduce the following symbols: 

(M \ I M+j\ M M M+j'^-^ M+j'^^^ 

E n Eh E ••• E E ••• E ■ 

^min ^max 

(5.17) 
Recall that a{j) are defined in (5.8). We define f3{j) by 

/?(j)=i-25[(j-l)/2s] (l<j<M). (5.18) 

For ^, fc == 1, 2, . . . , 2sm, we define the {£, k) element of M^'^'"'\{bj)2sm) by 

(m(2-)((5,.)2.™))^^ 

^f -<56,-M,fe {bi>M) 

\<5/3(fc,),/3(fe) -pCAfe, -wf"' +?7/2)/(A^sPtot(Aia(fc,))) {bi<M) 

Here, continuous variable /i, which is the argument of density ptot(M)i is 
evaluated at /ia(6f)i one of the "string centers" ^a of 2s-strings (5.7). 

We can rigorously derive a concise expression of correlation functions of 
the spin-s XXZ spin chain in the massless region: < C < 7r/2s for a large 
finite chain. Introducing (/3(A) — sinhA we have 

/ M \ (m+j\ 

Fi?,^+\{^,,3b})^C{{^u,Ju}) I] E IT E detM(2-)((6,)2,„) 

. n,e«- (ncl ^(Ac, - -f + ^) nr=7+i ^(^^^ - t^f^) 

x(-l) ^ 



ni<fe<£<25m "P^^hi - K + V) 



^ ' ^ '- + 0{1/Ns) 



ni<fe<£<2.m'/'K^"' -"^ 



(2«) _„,(2s)^, 

(5.20) 
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We remark that we derive (5.20) sending e to zero. Before taking the limit, 
inhomogeneous parameters WjS are generic due to small parameter e, and 
the sums over variables Cj in (5.20) are restricted up to M for all j. 

5.5. Multiple-integral representations of spin-s XXZ 
correlation function for arbitrary matrix elements 

In the thermodynamic limit: A^s — > oo, rapidities A;,^ with bi defined in 
(5.16), correspond to integral variables \i for i ~ 1,2, . . . ,2sm. For 1 < 
bt < M they are given by the Bcthc roots of 2s-strings (5.7), while for 
bi> M they are given by complete 2s-strings w^, " defined by (4.1). 

We define a(Aj) by a{\j) = 7 for an integer 7 with 1 < 7 < 2s, if A^ is 
related to integral variable fXj by Aj = fij — (7— l/2)?7, or if \j takes a value 
close to w^ " with /?(fc) = 7, where wj, ^ are part of complete 2s-strings 
(4.1). Here, variables Hj correspond to "string centers" of variables Aj. 

We define the (j, k) element of matrix S = S I {\j)2sm', {Wj ^ )2sm ) by 

Sj.k=p{\,-wf'^ +7l/2)5{a{X,),l3{k)), for j, fc = 1, 2, . . . , 2sm . 

(5.21) 
Here 5{a,(i) denotes the Kronccker delta, and wc recall (5.18) for /3(fc). 

(2s) 

Let Tj be a small contour rotating counterclockwise around A = w, . 
Since detS* has simple poles at A = w- ^ with residue l/27ri, we have 

00+ie poo — ie p 

detS{{Xk)2sni)dXi ^ detS{{Xk)2sni)dXi-(ti detS{{Xk)2sm)dXi 

-oo-\-ie J —oo—ie JT\ 

(5.22) 
For sets ol and ol^ with relation (5.16), wc define integral variables Xj for 
j G OL and A for j G ol^ , respectively, by the following: 

'-■^Jmax ' • • • ' '^Z . ' •^i".i" ' ^Jmax ) = ("^1 J ' • ' J Msm) ■ (5.23) 

Thus, from expression (5.20) of the correlation function in terms of a finite 
sum, we derive the multiple-integral representation as follows. 

F^n'^\{^k,]k})^ C{{i,,3k))y. 

dAi • • • / +•••+/ _ 1 dX^, 

-oo — -it^s+ie/ \J — oo-{-ie 




dXr- ■ 



-oo — i(^s—ie J \'^ —CO- 

xQ{{ej,e^}; Ai, . . . , X2sm) detS'(Ai, . . . , X2sm) 



oa—iC,s+itJ 




co — iQs—i^ \ 






dX2sm 


co-iCs-ieJ 






(5.24) 
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Here C^ = (2s - 1)C, r = r' + l, and Q{{€j,e'^}; Ai, . . . , \2sm) is given by 



ni<fc<£<2sm '/'(A£ - Afe + ry + e^^fc) 

a-ec. (nu ^(a; - -f ^ - ^) nr=7+i ^(^; - 4'^^) 



(5.25) 



ni<fc<£<2sm ¥'(w'f '■* - ^f '^) 

In the denominator we set ek.e — *e for Im(Xk — Xi) > and ek^i = — ie for 
Im{Xk — Xe) < 0, where e is infinitesinially small: |e| <C 1. Here, Im{a+ib) = 
b for real numbers a and 5. Here, for a , we recall (5.14). 

We evaluate a(Xj) in (5.24), replacing paths (— oo — i(j — 1)( ± ze, oo — 
z(7 - 1)C ± le) by (-(^ - i(7 - 1/2)C, (^ - z(7 - 1/2)C) for 7 = 1, 2, ... , 2s, 
respectively. The integrals over Xj for j > f do not change when e -^ C,/2. 

Thus, correlation functions (5.9) are expressed in the form of a single 
term of multiple integrals (5.24). 

We can derive the symmetric expression for the multiple-integral repre- 
sentations of the spin-s correlation function -Fm* {{ik,jk}) as follows.^ 

F!^'+\{^k,Jk})^C{{^,,J,})x 

1 T-r sinh2^(7r(6-6)/C) 



n 



ni<.<^<2. sinh™(/3 - a)rj ^JJ^^^ n -li UlU si^M^ - ^i + {r - j)v) 



y: (sgna)n / +•••+/ 



/oo — i(2s — l)(^-\-ie \ 



2sm 
X 



/oo— i(2s— l)t^— ie \ 
-oo-J(2s-l)C-«/ 



, . -' — 00 — ze 

j=r +1 ^ 



X (5({ej,ej}; Acti,. . ., Ao.(2sm))) ]_{ 



nl^i n^li^ sinh(Aj - Cb + /3?7) 



j^i nl^i cosh(7r(/^j - Cb)/C) 

^^— ^ n n sinh(7r(/^2s(a-l)+7 - M2.(6-l)+7)/C) • 



(5.26) 



Here Xj are given by A^ = /i^ — (/3(j) — 1/2)77 for j = 1, . . . , 2sm. 

It is straightforward to take the homogeneous limit: ^fc -^ 0. Here (sgn cr) 
denotes the sign of permutation <t € S2sm / {Smf''^ ■ 
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6. Derivation of finite-sum expression of spin-s XXZ 
correlation functions with arbitrary entries 

6.1. Fundamental commutation relations 

We now discuss briefly the derivation of (5.20), wliicli expresses tlie spin-s 
XXZ correlation functions with arbitrary entries in terms of the product of 
finite sums over the Bethe roots. 

Let Eat be the set of integers 1,2, ... ,N, i.e. Eat = {1, 2, . . . , N}. Recall 
definition (5.14) of a^ and that of integers Cj and c,. For a given set of 
Cj,Ca, we introduce Aj and A, by 

Aj = {b;l <b< M + 2sm, fe 7^ c^, c^ for fc < j} , 

a] ^ {b;l<b<M + 2sm, fe 7^ Cfe for fc < j, 6 7^ c), for fc < j}. (6.1) 

We define sets aj' and c{a^) as follows. 

aJ = {k; k < j, k e a^} , q:+ == {fc; fc < j, k G «+} , (6.2) 

c{aj) = {cfc; k e aj} , c{a+) = {c^; k e a+} . (6.3) 

We have 

Aj = TjM+2sm \ (c(a7) U c(a+)) , A^ == EM+2sm \ (c(q;J_|_i) U c{a.+)) . 

Let us denote by t the number of Cj {j G ol^) and c (7 G «+) such 
that Cj, c < M , for a given set of Cj and c . We express (5.17) as follows. 



(6.4) 



Here the sum over {cj,c,}t denotes the sums over Cj and c such that the 
number of c < M is fixed hy t — r. 

Suppose that Aq, for a = 1, 2, . . . , M give a set of solutions of the Bethe 
ansatz equations in the spin-1/2 case with Wj — w, *' '^ for j = 1, 2, . . . , L.^ 
Here Wj are inhomogeneous parameters. We set rapidities \m+j by 

Am+j — Wj , j ~ 1,2,. . . , 2sm . (6.5) 

We can show the fundamental commutation relations as follows. ^^ 
/ M \ 



2sm 


= 11 


M \ 

>; 


11 


/ M+j 

>; 


*=Mc,x;}, 


jea- 


\c,=l) 


jecx.+ 


\c,=l 



\{c{K)]t.{\m+,)---t^^_>J\ 



^Q = l 



M+2siyi) 



2sni 
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where d{X; {wf ^'"''11) and G^^^^^'y{{Xa)M+2sm) are given by 

k=l 
G{c„c'^}i^ir-- ,^M+2s,n)= n y.M+j ., « . 



i6"- V '■h^iMeA 






(6.6) 



6.2. Finite-sum expression of correlation functions for a 
finite chain 

We introduce disjoint subsets of a"'", a'^ and a~^, as follows. 

"j = {j; J e a+, 1 < c;. < M} , a+ = {j; J e a+, c;. > M} . (6.7) 

We define sets c{cy.^), c{a.y] and c{a.'^) as follows. 

c{ct^) = {ck; k e a"}, c(q:]") = {cfe; k e a]"}, c(q:;^) = {cfe; A: G a;^} . 

We define a sequence (6fc)t by a subsequence of {hk)2sm such that bk < M 
for fc = 1,2, ... ,t. We denote sequence {bk)2sm and (5fe)t as sets by b and 
bt, respectively, i.e. b ^ {&i, &2, • • • , fe2sm} and bt = {61, . . . , bt}. Here we 
note 6t = c(q:^) U c(q:j). We define sequence {bi.)2sm-t by a subsequence 
of {bk)2sm such that 6^, > M for fc = 1, 2, . . . , 2sm ~ t. We denote it as a 
set by 62sm-f Here we note b^^^_^ = c(a^).^ 

We define sets Z and iiT by Z = Sm \ bt and X = S2sm \ ^2sra-u 
respectively. We define a sequence (z(Q;))M-t by putting the elements of Z 
in increasing order: z(l) < z(2) < • • • < z(M — i) where Z = {z{i);i = 
1,2, . . . , M — t}, and a sequence (kj)^ by putting the elements of K in 
increasing order: ki < K2 < ■ ■ ■ < Kt where K = {kj; j — 1,2, ... , t}. 

We derive the spin-s correlation functions from those of the spin- 1/2 
case sending e to zero: 

J^i'^+HK Jj; (t^f "+')l) = C({z,, jU) liniFi^+'({6,,e;.}; (u;f ^^))i) . 

(6.8) 
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Applying (6.6) to (5.9) (or (5.13)) we have 

2s7n 

, (oi nfji* g(A.(.)) n;^i c{w,.)nf=-i Bjx^wmu b(vj\^) 

X02.™(iA.|M) ^^1 ^^_, ^(,^^^^) j^* ^^ ^(^^) ^M-. ^(,^^^^) j^. ^^ ^(,^)„^ 

2sm M 2sm ,, , /ttM+j-I /> ^ \, i n\ 

E\r^ TT TT y'(Aa - Wj) -p-r ( 1 lb=l,6eAj '/'I'^Cj ^b + Tj ) 

2^ 11 11 roQ, _7„.. ^'n\ 11 



t=r 



^ TT TT y(A2:(«) - Ag;) -|-|. -|-|- y(Aa - W^, + ?7) 
X dot I ($')-! (A,(„))M-t#(A,~)t) X 

X*' {{K(»))M-t#{w^,)u {K{»))M^t#{Xg^)t) 1 .(6.9) 



rM 



Here, 0„i({Aq}) == HjLi 110=1 ^(Aa - Wj), and matrix elements (* )„;, for 
a,b ~ 1,2, . . . , M are given by 



* ((A2(a))M-t#(Ar)t, (A;,(a))M-i#(wKjt;(Wfc)L), 

' a, 6 



*a,6((^-("))^^-*#(^bl)* 



for & < M - i 



in(ri) , (6.10) 

^^ ''^ for6 = A: + M - i (1 < fc < i) 



¥'(Az(a) - M^fcJ<^(A^(a) - w„, + ry) 
The matrix elements of the Gaudin matrix are given as follows 

*1,6 \{K(c))M-t^iX^)t\ (wfc)LJ = *l(a),z(6)((Aa)M; («^fe)L) 



</?(Aa - Ab + r\)if[\a - Ab - ry) V "^ "^(Aa - u'p)<p(Aa - Wp + ry) 

^^ </j(Aa - A^ + 77)(^(Aa -X^-v))' 
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For any positive integer Ng we can rigorously calculate (6.9) as follows. ^^ 
Proposition 6.1. 

F^i:\{e,.e^};{wf-'\)=Y^Y: n (-1) 
x(-ir"-* n n (-!)• n (-i) 

jeOL+, WeQ;+;£>j KeK;K+M<c'. 



X {-ly 



: Tl.ec- nU ^(Ac, - wt' '' + V) nr=7+i ^(Ac, - «.f ^ '^) 



Ill<k<e<2sm ^i^be - Afc^ + 7?) 
i.ll<k<l<2sm^\^k '^l I 

(6.12) 
We define matrix elements {j, k) of 0^/" ((&£)2sm) (1 < j < 2sm):^^ 
If 6, > M, f^if "\(6£)2.m)) - -<5b,-M,fc for fc = 1, 2, . . . , 2sm, 

\ / j,k 

if bj < M, there is an integer i such that bj = &i 

(</>!^^'"^((6£)2.™)) = ($')-i*'((A,(„))M-t#(A,-^)*, 

(Az(Q))M-t#(CK£)t)i+M-t,fc+M-t , for A: = 1,2, ... ,i, 
and (pM iibi)2srn)j^f,\^0 forfc == l,2,...,2sm-i. (6.13) 

We can show the following proposition.^^ 
Proposition 6.2. 

det (($y^*'((A^(„))M-t#(w«:Jt, (A^(„))M-t#(AfJt)) 

=dct0(^^™)((6,)2.™)(-i)^— * n (-1) 
X n n (-1) • n (-i) 1 • (6-i4) 

jeo;^ \eeoi.^;i>j k.£K;k,+m<c'. 
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When Ns is large enough, solving the integral equations for 
<^M " ((^'?)2sm): we Can show 

det^%'"'\{be)2srn) - det Af(2-™)((6,)2.™) + 0(l/iV.). (6.15) 

We thus obtain the finite-size spin-s XXZ correlation functions with arbi- 
trary entries (5.20). 
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